The unitary gas in an isotropic harmonic trap: symmetry properties and applications 
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We consider TV atoms trapped in an isotropic harmonic potential, with s-wave interactions of 
infinite scattering length. In the zero-range limit, we obtain several exact analytical results: mapping 
between the trapped problem and the free-space zero-energy problem, separability in hyperspherical 
coordinates, 50(2, 1) hidden symmetry, and relations between the moments of the trapping potential 
energy and the moments of the total energy. 
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Strongly interacting degenerate Fermi gases with two 
spin components are studied in present experiments with 
ultra-cold atoms p]]: by tuning the interaction strength 
between the atoms of different spin states via a Feshbach 
resonance, one can even reach the so-called unitary limit 
where the interaction strength in the s-wave chan- 
nel reaches the maximal amplitude allowed by quantum 
mechanics in a gas. More precisely, this means that the s- 
wave scattering amplitude between two particles reaches 
the value 



fk 



1 

ik 



(1) 



for the relative momenta k that are relevant in the gas, 
in particular for k of the order of the Fermi momentum 
kp of the particles. This implies that the s-wave scatter- 
ing length a is set to infinity (which is done in practice 
by tuning an external magnetic field). This also implies 
that k\r e \ <C 1, where r e is the effective range of the in- 
teraction potential, a condition well satisfied in present 
experiments on broad Feshbach resonances. 

The maximally-interacting gas defined by these con- 
ditions is called the unitary quantum gas Q|. It has 
universal properties since all the details of the interac- 
tion have dropped out of the problem. Theoretically, for 
spin 1/2 fermions with equal populations in the two spin 
states, equilibrium properties have been calculated in the 
thcrmodynamical limit in the spatially homogeneous case 
using Monte-Carlo methods; at finite temperature 
and at zero temperature with a fixed node approximation 
In practice, the unitary gases produced experimen- 
tally are stored in essentially harmonic traps, which rises 
the question of the effect of such an external potential. 
In this paper, we consider a specific aspect of this ques- 
tion: restricting to perfectly isotropic harmonic traps, 
but with no constraint on the relative spin populations, 
we show that the unitary quantum gas admits interesting 
symmetry properties that have measurable consequences 
on its spectrum and on the many-body wavefunctions. 
These properties imply that there is a mapping between 
the A^-body eigenfunctions in a trap and the zero-energy 
iV-body eigenfunctions in free space; the iV-body prob- 
lem is separable in hyperspherical coordinates; and there 
exist relations between the moments of the trapping po- 



tential energy and those of the total energy at thermal 
equilibrium. 

A unitary Bose gas was not produced yet. This is re- 
lated to the Efimov effect 6] : when three bosons interact 
with a short range potential of infinite scattering length, 
an effective three-body attraction takes place, leading in 
free space to the existence of weakly bound trimers. This 
effective attraction generates high values of k so that the 
unitarity condition Eq. Q is violated. It also gives a short 
lifetime to the gas by activating three-body losses due to 
the formation of deeply bound molecules [3, l& li| ■ In an 
isotropic harmonic trap, for three bosons, there exist efi- 
movian states 0,^1, but there also exist eigenstates not 
experiencing the Efimov effect |lCl [l2j . These last states 
are universal (in the sense that they depend only on h, 
the mass m and the trapping frequency u) and they are 
predicted to be long-lived [13 • The results of the present 
paper apply to all universal states, fermionic or bosonic, 
but do not apply to the efimovian states. 



I. OUR MODEL FOR THE UNITARY GAS 

The physical system considered in this paper is a set 
of A^ particles of equal mass m (an extension to dif- 
ferent masses is given in Appendix [SJ . The particles 
are of arbitrary spin and follow arbitrary statistics; the 
Hamiltonian is supposed to be spin-independent so that 
the A^-body wavefunction ip that we shall consider corre- 
sponds to a given spin configuration 0. The particles 
are trapped by the same isotropic harmonic potential of 
frequency u. We collect all the positions ?i of the parti- 
cles in a single 3 A" component vector: 



X = (fx, . . .,f N ). 



Its norm 



X = \\X\ 



N 



(2) 



(3) 



is called the hyperradius. We will also use the unit vector 



X/X 



(4) 
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(which may be parametrized by 3jV— 1 hyperangles). The 
coordinates (X, ft) are called hyperspherical coordinates 
[l4j . The total trapping potential energy simply writes 



trap 



(5) 



The interaction between the particles is assumed to 
be at the unitary limit defined in Eq.Q; one can then 
replace the interaction by contact conditions on the N- 
body wavefunction (see e. g. 0, 0] and references 
therein): when the distance — \\fj —fi\\ between par- 
ticles i and j tends to zero, there exists a function A such 
that 



A(Rij,{f k :k^i,j}) 



0{n 



(6) 



where Rij = (fi + fj)/2 is the fixed center-of-mass posi- 
tion of particles i and j, and {fk : k ^ are the po- 
sitions of the other particles. When none of the particle 
positions coincide, the stationary wavefunction if> solves 
Schrodinger's equation, Hip = Etp, with the Hamiltonian 



H 



2m x 



1 



mu 2 X 2 . 



(7) 



This model is expected to be exact for universal states in 
the limit of a zero range of the interaction potential . 



II. SCALING PROPERTIES OF THE TRAPPED 
UNITARY GAS 

A. What is scale invariance ? 

A fundamental property of the contact conditions 
Eq.@) is their invariance by a rescaling of the spatial 
coordinates. More precisely, we define a rescaled wave- 
function tpx by 



Mx) = ^(x/\) 



(8) 



where A > is the scaling factor. Then, if -0 obeys 
the contact conditions, so does ip\ for any A. Note that 
this property holds only because the scattering length is 
infinite (for a finite value of a, 1/rij in Eq.JBJ) would be 
replaced by — 1/a, which breaks scale invariance). 

In free space (that is for u> = 0), the scale invariance 
implies the following property: if ip is an eigenstate of 
energy E, then ipx is an eigenstate of energy E/X 2 for any 
A [lTj ■ This implies the absence of bound states in free 
space: otherwise the scaling transform would generate 
a continuum of states which are square integrable (after 
elimination of the center of mass variables), and this is 
forbidden for a Hermitian problem |isj|. 

When E = 0, one finds (see Appendix |B|) that the free 
space eigenstates can be assumed to be scale-invariant, 
i.e. there exists an exponent v such that 



(9) 



Taking the derivative of this relation with respect to A in 
A = 1, this shows that ip is an eigenstate of the dilatation 
operator, 



D = X-d^ 



(10) 



with the eigenvalue v. This result is interesting for sec- 
tion EH 

The presence of a harmonic trap introduces the har- 
monic oscillator length scale a^o = ijh/mw, so that the 
eigenstates cannot be scale- invariant as in Eq.@. How- 
ever, if ip obeys the contact condition, so do the tpx's: as 
we shall see, this allows to identify general properties of 
the eigenstates in the trap. 



B. Scaling solution in a time dependent trap 

We now assume that the trap frequency to, while keep- 
ing a fixed value w(0) for all times t < 0, has an arbitrary 
time dependence at positive times. Let us assume that, 
at t < 0, the system is in a stationary state of energy 
E. Then at positive times the wavefunction of the sys- 
tem will be deduced from the t = wavefunction by the 
combination of gauge and scaling transform |l9j |: 

-iEr(t)/h 

i,(X,t) = e imX W w ^(X/A(i),0) (11) 

where the time dependent scaling parameter obeys the 
Newton-like equation 



A 



W 2 (0) , 



A 3 



(12) 



with the initial conditions A(0) = 1,A(0) = 0. We also 
introduced an effective time r given by 



r(t) 



* dt' 



o A 2 (t')' 



(13) 



This result may be extended to an arbitrary initial 
state as follows: 



rP(X,t) 



1 



A(i) 3W / 2 



^ k ^l 2hx ^i>{X/\{t),r{t)). 



(14) 

where ip evolves with the t < Hamiltonian (i.e. in the 
unperturbed trap of frequency w(0)). 

As shown by Rosch and Pitaevskii |2Cl |. the existence 
of such a scaling and gauge time dependent solution is 
related to a SO(2,l) hidden symmetry of the problem. 
This we rederive in the next subsection. 



C. Raising and lowering operators, and SO(2,l) 
hidden symmetry 

We consider the following gedanken experiment: one 
perturbs the gas in an infinitesimal way by modifying the 



3 



trap frequency in a time interval < t < tf. After the 
excitation period (t > tf), the trap frequency assumes its 
initial value oj(0). The scaling parameter then slightly 
deviates from unity, X(t) = 1 + SX(t) with \SX\ <C 1. 
Linearizing the equation of motion Ea. (|12l) in <5A, one 
finds that <5A oscillates as 



SX(t) 



(15) 



where we set to = uj(0) to simplify the notation. The 
gedanken experiment has therefore excited an undamped 
breathing mode of frequency 2ui (2£j . 

We now interpret this undamped oscillation in terms 
of a property of the ./V-body spectrum of the system. 
Expanding Ea. Hllfl to first order in 5X(t) leads to 



+ e * e -HE-2Hu)t/h L ] o) + 0( e 2 ) (16) 



(the phase a depends on the details of the excitation pro- 
cedure). This reveals that the initial stationary state E 
was coupled by the excitation procedure to other station- 
ary states of energies E ± 2hui. Remarkably, the wave- 
function of these other states can be obtained from the 
initial one by the action of raising and lowering operators: 

3iV H 
L+ = +—+b+ — -m,ujX 2 /h (17) 
Z Tioj 

L_ = -— - D + — - ttiluX 2 /%. (18) 
2 nuo 

Repeated action of L + and L_ will thus generate a ladder 
of eigenstates with regular energy spacing 2huj. 

The hidden SO(2,l) symmetry of the problem then re- 
sults from the fact that H , L + and L_ have commutation 
relations equal (up to numerical factors) to the ones of 
the Lie algebra of the SO(2,l) group, as was checked in 
M: 



[H,L, 
[H,L. 

[L + ,L_] = -4 



= 2fkoL + 
= -27ui)L- 
H 
fi/jj 



(19) 
(20) 

(21) 



From the general theory of Lie algebras, one may form 
the so-called Casimir operator which commutes with all 
the elements of the algebra, that is with H and L±; it is 
given by [20|: 



C = H z 



! (L+L- +L-L-I 



Consider a ladder of eigenstates; as we will show later, 
the hermiticity of H implies that this ladder has a ground 
energy step, of value Eq. Within this ladder, the Casimir 
invariant assumes a constant value, 



D. Virial theorem 

Another application of the existence of raising and 
lowering operators is the virial theorem for the unitary 
gas. For a given eigenstate of H of energy E and real 
wavefunction ip, L^\^p) is either zero (if ip is the ground 
step of a ladder) or an eigenstate of H with a differ- 
ent energy. Assuming that H is hermitian, this implies 
(ip\L-\ip) = 0, and leads to H3: 

<VW)=2(i/>|ff trap k/>)- (24) 
At thermodynamical equilibrium, one thus has 



(H) = 2 (H tIap ) , 



(25) 



that is the total energy is twice the mean trapping po- 
tential energy. 

This virial theorem is actually also valid for an 
anisotropic harmonic trap (this result is due to Frederic 
Chevy). One uses the Ritz theorem, stating that an 
eigenstate of a hermitian Hamiltonian is a stationary 
point of the mean energy. As a consequence, the function 
of A 



E(X) 



(V'aIV'a) 

= A" 2 (i/>\H - H tiap \ip) + A 2 MfltfpM (26) 



satisfies (dE /d\)(\ = 1) = 0, which leads to the virial 
theorem. This relics simply on the scaling properties of 
the harmonic potential, irrespective of its isotropy. 

The proportionality between (H) and {H tla , p ) resultin 
from the virial theorem was checked experimentally 



.ting 



III. MAPPING TO ZERO-ENERGY 
FREE-SPACE EIGENSTATES 

Usually, the presence of a harmonic trap in the exper- 
iment makes the theoretical analysis more difficult than 
in homogenenous systems. Here we show that, remark- 
ably, the case of an isotropic trap for the unitary gas can 
be mapped exactly to the zero-energy free-space prob- 
lem (which remains, of course, an unsolved many-body 
problem). 

More precisely, all the universal TV-body eigenstates 
can be put in the unnormalized form: 



|Vv, g > = (£+)* e -* 2/2a 'H€> 



(22) an d have an energy 



EL 



- 2q + 2,N/2)fiuj 



(27) 



(28) 



where q is a non-negative integer, L+ is the raising opera- 
tor defined in Eq. I|17(l , and ip® is a zero-energy eigenstate 
of the free-space problem which is scale-invariant: 



C = E (E -2hu). 



(23) 



(29) 
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for all real scaling parameter A, v being the real scaling 
exponent 

We also show that the reciprocal is true, that is each 
zero-energy free-space eigenstate which is scale-invariant 
with a real exponent v generates a semi-infinite ladder of 
eigenstates in the trap, according to Ea. H27l28(l . 

We note that Eq.(l28|l generalizes to excited states a 
relation obtained in [24| for the many-body ground state. 



A. From a trap eigenstate to a free-space 
eigenstate 

We start with an arbitrary eigenstate in the trap. By 
repeated action of L_ on this eigenstate, we produce a 
sequence of eigenstates of decreasing energies. According 
to the virial theorem Eq. I|24(l , the total energy of a univer- 
sal state is positive, since the trapping potential energy is 
positive. This means that the sequence produced above 
terminates. We call ip the last non-zero wavefunction 
of the sequence, an eigenstate of H with energy E that 
satisfies L-\ip) = 0. To integrate this equation, we use 
the hyperspherical coordinates (X, n) defined in Eq. H3I4J) . 
Noting that the dilatation operator is simply D — X ■ dx 
in hyperspherical coordinates, we obtain: 

ip(X ) - e - x2 / 2 < x E ^- 3N ^ 2 f(n). (30) 



Then one defines 



V>°(A) 



(31) 



One checks that this wavefunction obeys the contact con- 
ditions Eq.jnj, since X 2 varies quadratically with at 
fixed Rij and {fk,k ^ '0° is then found to be 

a zero-energy eigenstate in free space, by direct inser- 
tion into Schrodinger's equation. But one has also from 
Ea. CTTll : 



(32) 



so that ip° is scale- invariant, with a real exponent v re- 
lated to the energy E by Ea. H28|) . This demonstrates 
Ea. (|27l28|l for q = 0, that is for the ground step of each 
ladder. 

One just has to apply a repeated action of the raising 
operator L+ on the ground step wavefunction to generate 
a semi-infinite ladder of eigenstates: this corresponds to 
q > in Ea. l|27l28|l . Note that the repeated action of 
L + cannot terminate since L + \ip) = for a non-zero tp 
implies that ip is not square-integrable. 



resulting wavefunction is an eigenstate of the Hamilto- 
nian of the trapped system, obeying the contact condi- 
tions [25|. Applying L + then generates the other trap 
eigenstates. 



C. Separability in hyperspherical coordinates 

Let us reformulate the previous mapping using the hy- 
perspherical coordinates ( X, n) defined in Eq. l|3l4|l . A 
free-space scale- invariant zero-energy eigenstate takes the 
form ip°(X ) = X v fv(n), and the universal eigenstates in 
the trap have an unnormalized wavefunction 

i> v , q (X) = X»e- x2 / 2 <L^ N / 2 Hx 2 /a 2 )U(n), 

(33) 

where L q is the generalized Laguerre polynomial of de- 
gree q. This is obtained from the repeated action of L + 
in Eq. (|27|l and from the recurrence relation obeyed by 
the Laguerre polynomials: 



(q+l)L^(u)-(2q+S+l-u)L q s \u) + (q+s)L q s l l (u) = 0. 

(34) 

We have thus separated out the hyperradius X and the 
hyperangles n. The hyperangular wavefunctions 
and the exponents v are not known for N > 4. However, 
we have obtained the hyperradial wavefunctions, i.e. the 
X dependent part of the many-body wavefunction. A 
more refined version of these separability results can be 
obtained by first separating out the center of mass (see 
Appendix Q, but this is not useful for the next Section. 



-(*) 



IV. MOMENTS OF THE TRAPPING 
POTENTIAL ENERGY 

A. Exact relations 

As an application of the above results, we now obtain 
the following exact relations on the statistical properties 
of the trapping potential energy, relating its moments to 
the moments of the full energy, when the gas is at thermal 
equilibrium. For the definition of the trapping potential 
energy, see Eq.©. 

At zero temperature, its moments as a function of the 
ground state energy Eq are given by: 

((fftr ap ) n ) = E Q (E + huo)... (E + (n - l)hu) IT. (35) 

At finite temperature T, the first moment is given by the 
virial theorem 



B. From a free-space eigenstate to a trap eigenstate 

The reciprocal of the previous subsection is also true: 
starting from an arbitrary zero-energy free-space eigen- 
state that is scale-invariant, one multiplies it by the 
Gaussian factor exp(— X 2 /2a 2 10 ), and one checks that the 



(fftrap) - (H)/2 

and the second moment by 



((-Htrap) ) - 



(H 2 ) + {H)fiuj ■ cotanh 



(36) 



/4. 



(37) 
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Derivation 



CONCLUSION 



The zero temperature result Eg. (jr{5|) follows directly 
from Eg. H33|) : for q = 0, the Laguerre polynomial is con- 
stant so that the probability distribution of X is a power 
law times a Gaussian; the moments are then given by in- 
tegrals that can be expressed in terms of the T function. 

For finite T, the idea of our derivation is the following: 
the hyperradial part of the ./V-body wavefunction ip u _ q is 
known from Eq. Ij33|l ■ and thus the probability distribu- 
tion of X in the state \^ v , q ) is known, in terms of is, q. 
While the thermal distribution of q is simple, the one of 
v is not, but v is related to the total energy by Ea. (|28|l . 

We will need the intermediate quantities: 



J-due-^L q %(u)L q s >(u) 



-00, 



B 



f™du> 



L q s) iu) 



(38) 



where s > 0; n, q are non-negative integers; and p is an 
integer of arbitrary sign. These quantities can be calcu- 
lated with the n — 'initial' condition B$ p — <5o, P and 
the recurrence relation 



In this paper we have derived several exact properties 
of the unitary gas in an isotropic harmonic trap. The 
spectrum is formed of ladders; the steps of a ladder are 
spaced by an energy 2fuv, and linked by raising and low- 
ering operators. This allows to map the trapped prob- 
lem to the free-space one. The problem is separable in 
hyperspherical coordinates. This allows to derive exact 
relations between the moments of the trapping potential 
energy and the moments of the total energy. The rela- 
tion between the first moments is the virial theorem; the 
relation between the second moments may be useful for 
thermometry, as will be studied elsewhere. 

We thank F. Chevy, J. Dalibard, S. Nascimbene, L. 
Pricoupenko, J. Thomas and C.-F. Vergu for very useful 
discussions. LKB is a Unite de Recherche de I'ENS et 
de VUniversite Paris 6, associee au CNRS. Our research 
group is a member of IFRAF. 



APPENDIX A: EXTENSION TO PARTICLES 
WITH DIFFERENT MASSES 



B n+ i, p = -(q+p+ l)B n , p+1 + [2{q+p) + s + l]B ntP 
-(q + p + s)B n ^ 1 (39) 

which follows from the recurrence relation Eq. I|34|) on La- 
guerre polynomials. 

This allows to calculate the moments of the trapping 
energy in the step q of a ladder of exposant v, using 
Eq.(E2J|: 



(^, q \X 2n \j>„, q ) 

{^v,q\^u,q) 



&n,0 a ho . 



Here we have set 



s = v-\ + 3N/2 



(40) 



(41) 



in accordance with Ea. l|33|l . 

Assuming thermal equilibrium in the canonical ensem- 
ble, the thermal average can be performed over the sta- 
tistically independent variables q and s. The moments of 
q are easy to calculate, because of the ladder structure 
with equidistant steps: 



(<f ) = 



+°° n n -qhul/k B T 



E + oo 



E + oo 
<?=0 1 



-qhuj/kBT 



(42) 



The moments of s are not known exactly but they can 
be eliminated in terms of the moments of the total en- 
ergy E and of the moments of q using the relation 
E = (s + 1 + 2q)hui. This leads to the exact relations 
(|36I37[1 . This method in principle allows to calculate re- 
lations for moments of arbitrary given order, but the al- 
gebra becomes cumbersome. 



All our results remain valid if the particles have dif- 
ferent masses mi,...,mjv; provided that the trapping 
frequency to remains the same for all the particles. We 
define a mean mass: 



mi 



m N 



N 



(Al) 



The definition of X and X, given by Eq. l|2l3|l for equal 
masses, has to be generalized to: 



X 



X 



n, ■ 



m N 



-r N 



\X\ 



N 



1 2 
-77 . 



(A2) 



(A3) 



With this new definition of X, the trapping potential 
energy is still given by Eq. © . 

In the definition of the zero-range model, the con- 
tact conditions Eq.® remain unchanged, except that the 
fixed center of mass position of particles i and j is now 
Rij = {rriifi + rii£fj)/(mi + rrij). 

In Appendix |C] the center of mass position has to be 
redefined as 



C = 



m N r N ) 



(mi + . . . + tow) 
and the internal hyperangular coordinates become 



R 



n 




(A4) 

(A5) 
(A6) 
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With these modified definitions, all the results of this 
paper remain valid. 

APPENDIX B: SCALE INVARIANCE OF THE 
ZERO-ENERGY FREE-SPACE EIGENSTATES 

In this appendix, we show that the zero-energy free- 
space eigenstates of the Hamiltonian may be chosen as 
being scale- invariant, that is as eigenstates of the dilata- 
tion operator D, under conditions ensuring the hermitic- 
ity of the Hamiltonian. 

Consider the zero-energy eigensubspace of the free- 
space Hamiltonian. This subspace is stable under the ac- 
tion of D. If one assumes that D is diagonalizable within 
this subspace, the corresponding eigenvectors form a 
complete family of scale invariant zero-energy states. If 
D is not diagonalizable, we introduce the Jordan normal 
form of D. 

Let us start with the case of a Jordan normal form of 
dimension two, written as 

Mat(D)=^^, (Bl) 

in the sub-basis \e\), |e 2 ). The ket |ex) is an eigenstate of 
D with the eigenvalue v. We assume that the center of 
mass motion is at rest, with no loss of generality since it is 
separable in free space. Using the internal hyperspherical 
coordinates (R, ft ) defined in Appendix [U] we find that 
D reduces to the operator ROr. Integrating Rdaei = ve% 
leads to 

e x {X) =R' / cf>i(n). (B2) 

The ket e 2 ) is not an eigenstate of D but obeys Rdsfii = 
vc2 + &x, which, after integration, gives 

e 2 (X) =R v logR<f> 1 (ti) + R"Mti)- (B3) 

One can assume that <pi and </> 2 are orthogonal on the 
unit sphere (by redifining e 2 and 2 ). It remains to 
use the fact that both e\ and e 2 are zero-energy free- 
space eigenstates. From the form of the Laplacian in 
hyperspherical coordinates in d — 3N — 3 dimensions, 
see Ea. (JC5|l . the condition A^ei = leads to 

Tj^i = -u(u + d-2)4> 1 . (B4) 

The condition A^e 2 = then gives Txfo = —v(y + d — 
2)02 — (2^ + d — 2)0i, which leads to the constraint [2(| : 

v = 1 - d/2. (B5) 

At this stage, for this 'magic' value of v, it seems that 
there may exist non scale-invariant zero-energy eigen- 
states. 

To proceed further, one has to check for the hermiticity 
of the free space Hamiltonian. This requires a reasoning 



at arbitrary, non zero energy. We use the fact that the 
following wavefunction obeys the contact conditions, 

ip(X) — u(R)R v <j)i(n), (B6) 

where u(R) is a fonction with no singularity, except 
maybe in R = [^j. Using again the expression of 
the Laplacian in internal hyperspherical coordinates, one 
finds that ij) is an eigenstate of the free-space Hamiltonian 
if u(R) is an eigenstate of 

h = -^(d 2 R + R- 1 d R ). (B7) 

One checks that Hermiticity of the free-space Hamilto- 
nian for the wavefunction ip implies hermiticity of h for 
the 'wavefunction' u(R). Note that h is simply the free- 
space Hamiltonian for 2D isotropic wavefunctions. It is 
hermitian over the domain of wavefunctions u(R) with a 
non-infinite limit in R = 0. Including the ket |e 2 ) in the 
domain of the iV-body free-space Hamiltonian amounts 
to allowing for 'wavefunctions' u(R) that diverge as log R 
for R — > 0: this breaks the Hermiticity of h, since this 
leads to a (negative energy) continuum of square inte- 
grable eigenstates of h, 

u K (R) = K (kR) (B8) 

with eigenenergy —fh 2 K 2 /2m, for all k > 0. Here Ko(x) is 
a modified Bessel function of the second kind. Hermitic- 
ity may be restored by a filtering of this contiunuum |29j | , 
adding the extra contact condition u(R) = log(i?//)+o(l) 
for R — ► 0, but the introduction of the fixed length I 
breaks the universality of the problem and is beyond the 
scope of this paper We thus exclude e 2 from the domain 
of the Hamiltonian. 

This discussion may be extended to Jordan forms of 
higher order. For example, a Jordan form of dimension 
3 generates a ket |es) such that (D — v)ez — e 2 . But e 2 
must be excluded from the domain of the Hamiltonian 
by the above reasoning. Since we want the domain to be 
stable under D, must be excluded as well. 

As a conclusion, to have a free-space A^-body Hamil- 
tonian that is both hermitian and universal (i.e. with 
a scale-invariant domain) forces to reject the non scale- 
invariant zero-energy eigenstates, of the form Eq.JBJJ. 

APPENDIX C: SEPARABILITY IN INTERNAL 
HYPERSPHERICAL COORDINATES 

We develop here a refined version of the separability 
introduced in subsection llll CI First, we separate out the 
center of mass coordinates. Then we obtain the separa- 
bility in hyperspherical coordinates relative to the inter- 
nal variables of the gas, which allows to derive an effective 
repulsive iV — 1 force and to get a lower bound on the 
energy slightly better than the one E > ensuing from 
the virial theorem. 
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Let us introduce the following set of coordinates: 



is the position of the center of mass (CM); 



(CI) 



R 



N 



(C2) 



is the internal hyperradius; and 

'rx-C rjv-C" 



fl = 



R 



R 



(C3) 



is a set of dimensionless internal coordinates that can be 
parametrized by 3iV — 4 internal hyperangles. In these 
coordinates, the Hamiltonian decouples as H = Hqm + 
i?int with 



Hcm 

Hint 



— — + l-Nmuj 2 C 2 
2Nm c 2 



2 m 
1 



3iV 



+ ^muj 2 R 2 



(C4) 



(C5) 



where X 1 ^ is the Laplacian on the unit sphere of dimen- 
sion 37V — 4. The contact conditions do not break the 
separability of the center of mass valid in a harmonic 
trap, so that the stationary state wavefunction may be 
taken of the form 



VPO = VCM(C) Vint i 



(C6) 



One can show 30] that there is separability in internal 
hyperspherical coordinates: 



v-intOM) = #(ij)0(n; 



(C7) 



This form may be injected into the internal Schrodinger 
equation 



-ff hit Vint = -EjntVin 



(C8) 



One finds that 0(ft ) is an eigenstate of with an eigen- 
value that we call —A. Note that the contact conditions 
Eq. © put a constraint on 0(fi* ) only j27j • The equation 
for $(i?) reads: 



3N 



R 



Or)® 



( h 2 A l 2 ,\ 
V 2^ + 2™ * 2 J $ 
= EM*. (C9) 



A useful transformation of this equation is obtained by 
the change of variable: 



$(R) = i? (5 - 3JV) / 2 F(i?), 



(CIO) 



resulting in 



^( 9 R+^r)f (^JL ± m ^lf)F 



f^s 2 , 1 



V2mi? 2 2 
= E int F{R), (CM) 



where sr is such that 



.si = A 



f3N 



V 2 



(C12) 



Formally, the equation for F is Schrodinger 's equation 
for a particle of zero angular momentum moving in 2D 
in a harmonic potential plus a potential oc s 2 R / R 2 . 

For s 2 R > 0, one can choose > 0. Assuming 
that there is no A^-body resonance, F(R) is bounded for 
R -> 31]. The eigenfunctions of Ea. (|CTll can then be 
expressed in terms of the generalized Laguerre polynomi- 
als: 



F(R) 



R SR L s q R [R 2 /al ]e- R /2a *° 



with the spectrum: 



E- m t = (sr + 1 + 2q)huj. 



(C13) 



(C14) 



This gives a lower bound on the energy of any universal 
A'-body eigenstate: 



E > -Klo 
~ 2 



(C15) 



for N > 2 and in the absence of a Af-body resonance. 

For a complex s R , the effective 2D Hamiltonian is not 
hermitian and this case has to be discarded. For s 2 R < 0, 
Whittaker functions are square integrable solutions of the 
effective 2D problem for all values E- m t so that, again, the 
problem is not hermitian. One may add extra boundary 
conditions to filter out an orthonormal discrete subset (as 
was done for A^ = 3 bosons (a, UJJ, liililJ] ) but this breaks 
the scaling invariance of the domain and generates non- 
universal states beyond the scope of the present paper. 

To make the link with the approach of Section llTTl we 
note that 



F(R) = R SR 



(C16) 



is a solution of the effective 2D problem (|C11(I for uj — 
0, £\nt = 0. Thus a solution of the internal problem 
Eci. (|C8(l at zero energy in free space is given by 

Vint(i?, ft) = i? (5 - 3JV)/2+S Mft). (C17) 

Multiplying this expression by C l Y l m (C/C), one recovers 
the VS's of Section Ulil with 

5 -37V 



+ Sr + I. 



(C18) 
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